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Leon  Clancy  Receives 
NASA  Awards 

Leon  M,  Clancy,  System  Manager,  was  awarded  a 
NASA  Public  Service  Medal  in  a  ceremony  held  at 
Langley  Research  Center  on  April  3,  1995.  This 
award  is  the  highest  honor  which  NASA  confers  to 
non-government  individuals,  and  was  given  For  sig¬ 
nificant  contributions  in  developing  a  low-cost  com¬ 
munications  network  that  allows  a  practical  High 
Performance  Computing  and  Communications  Ed¬ 
ucation  Outreach  Program  at  NASA  Langley  Re¬ 
search  Center.  Mr.  Clancy  is  also  a  member  of  the 
HPCCP  K-12  Educational  Outreach  Team  which 
received  a  NASA  Group  Achievement  Award  for 
its  work  in  developing  Internet-based  educational 
programs  for  school  systems  in  eastern  Virginia. 


Best  Mentorship  Award 

Shishir  Mehrotra  (pictured  on  the  right)  was  given 
the  Best  Mentorship  Project  award  by  the  New 
Horizons  Governor’s  School.  ICASE  Staff  Scientist 
David  Banks  (left)  served  as  his  mentor.  The  men¬ 
torship  program  involves  some  70  high  school  stu¬ 
dents  from  the  region.  Students  spend  at  least  five 
hours  a  week  with  a  mentor  at  a  science-oriented  in¬ 
stitution.  They  propose  projects,  conduct  research, 
and  then  present  their  results  at  the  end  of  the 
school  year. 

Shishir  studied  3D  computer  graphics  and  scien¬ 
tific  data  visualization.  He  wrote  software  to  ren¬ 
der  images  of  3D  objects,  and  he  tested  high-end 
visualization  systems  to  determine  their  usefulness 
for  analyzing  multivariate  and  unstructured  scien¬ 
tific  data.  _ 


Left  to  right:  A.  Thomas  Young,  Leon  M.  Clancy, 
and  Paul  F.  Holloway. 


Left  to  right:  David  Banks  and  Shishir  Mehrotra. 
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Use  of  the  World  Wide  WEB 
at  ICASE 

Leon  M.  Clancy* 

Introduction 

In  July  of  1993,  ICASE  began  to  use  the  World 
Wide  Web  (WWW).  The  Web  is  a  collection  of 
information  servers  on  the  Internet  offering  Hy¬ 
pertext  documents  written  in  Hypertext  Markup 
Language  (HTML).  The  Internet  includes  Com¬ 
mercial,  Governmental,  and  Educational  sites,  all 
of  which  may  participate  in  the  Web  by  making 
HTML  documents  available  through  a  Hypertext 
Transfer  Protocol  server  (HTTP).  Both  HTTP  and 
HTML  are  open  and  evolving  non-proprietary  stan¬ 
dards  offered  for  most  computers. 

Hypertext  documents  incorporate  a  structured 
method  for  defining  the  format  of  a  document  in 
a  system-independent  way.  The  document  may  in¬ 
clude  formatted  text,  graphics,  and  links  to  other 
Hypertext  documents.  Hypertext  documents  any¬ 
where  on  the  Web  are  referenced  using  a  Uni¬ 
form  Resource  Locator  (URL).  An  example  of  a 
URL  is  the  reference  to  the  ICASE  home  page: 
http:/ /www.icase.edu. 

The  large  number  of  information  sites  on  the 
Web  created  the  need  for  sophisticated  Web  brows¬ 
ing  and  searching  tools.  ICASE  initially  chose  Mo¬ 
saic,  one  of  the  earliest  graphical  Web  browsers 
available  for  all  of  our  computers.  Mosaic  is  a 
license-free  application  available  from  the  National 
Center  for  Supercomputing  Applications  (NCSA). 
It  can  be  accessed  on  the  WEB  at: 

http:/ /www. yahoo.com/Computers/World_Wide_Web/Browsers. 

Mosaic  offers  an  extremely  friendly  and  powerful 
interface  to  the  Web  which  is  easily  used  by  ICASE 
staff,  visitors,  and  consultants  to  obtain  up-to-date 
information  about  ICASE.  Since  the  scope  of  the 
Web  is  worldwide,  anyone  using  a  browser  is  also 
aUowed  to  access  ICASE  information  deemed  to  be 
of  interest  to  the  general  public. 

Designing  the  ICASE  Home  Page 

Web  browsing  usually  starts  with  an  anchor  point 
known  as  the  Home  Page.  The  ICASE  home  page 
was  created  to  make  available  the  who,  what,  when, 
where  and  why  of  ICASE,  in  addition  to  our  main 
research  output  of  ICASE  reports. 

*Leon  Clancy  is  System  Manager  at  ICASE.  He  joined 
ICASE  in  1989. 


The  ICASE  home  page  is  designed  in  such  a 
way  that  information  is  provided  and  maintained 
by  the  responsible  staff  members.  By  structuring 
the  home  page  to  integrate  the  efforts  of  the  staff, 
maintenance  of  the  home  page  is  divided  into  man¬ 
ageable  portions,  while  retaining  a  unified  appear¬ 
ance  and  direction.  This  avoids  the  need  for  a  staff 
member  dedicated  to  Web  maintenance.  Software 
tools,  often  Web  based,  are  provided  to  the  main- 
tainers  so  that  the  input  text  may  be  converted  to 
HTML  and  incorporated  into  the  home  page. 

Building  the  ICASE  Home  Page 

ICASE  technical  reports  form  the  largest  body  of 
information  we  make  available  on  the  Web.  Since 
all  reports  are  approved  for  unlimited  distribution, 
an  ICASE  report  can  be  distributed  via  the  Web 
to  a  world-wide  audience. 

In  cooperation  with  NASA/Langley,  ICASE  re¬ 
port  abstracts  are  full-text  indexed  and  integrated 
with  the  Langley  Technical  Report  Server  (LTRS). 
The  ICASE  Technical  Report  Server  is  used  to 
search  and  retrieve  reports  from  our  home  page. 
Our  part-time  technical  editor  is  responsible  for 
final  editing  and  reproduction  of  ICASE  reports. 
A  family  of  scripts  was  written  that  automate  the 
process  of  adding  the  final  report  to  the  ICASE 
Technical  Report  Server.  About  163  ICASE  re¬ 
ports  are  now  available  through  the  Web,  and 
about  7100  electronic  copies  have  been  downloaded 
through  our  server  since  December  of  1993. 

Another  interesting  use  of  the  web  at  ICASE 
is  the  Colloquium  announcement.  The  announce¬ 
ment  itself  has  a  simple  notice  consisting  of  the 
speaker’s  name,  title,  date,  time,  location,  and  ab¬ 
stract.  The  overall  operation  calls  for  announce¬ 
ment  generation  and  mailing  to  our  colloquium  list 
about  a  week  before  the  talk.  A  colloquium  is 
added  to  the  system  by  completing  a  web-based 
form  that  includes  this  information.  On  the  morn¬ 
ing  of  the  colloquium,  a  short  reminder  message 
is  then  mailed  to  a  list  of  people  who  have  asked 
to  receive  colloquium  announcements.  At  the  time 
the  announcement  is  generated,  it  is  converted  to 
HTML  and  placed  in  the  ICASE  Colloquium  sec¬ 
tion  of  the  home  page.  The  announcement  remains 
there  until  the  day  after  the  colloquium  has  been 
given,  at  which  point  it  is  moved  to  the  Past  Collo- 
quia  section.  There,  past  announcements  are  view¬ 
able  in  reverse  chronological  order. 
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There  are  a  wide  variety  of  web  browsers,  each 
of  which  has  slightly  different  features  and  capabil¬ 
ities.  One  design  goal  for  the  ICASE  home  page 
has  been  to  automatically  recognize  which  browser 
is  being  used  and  provide  a  version  of  the  ICASE 
home  page  that  is  adapted  to  the  requirements  of 
that  browser.  For  instance,  if  the  user  has  a  text- 
only  browser,  then  it  is  important  to  provide  suit¬ 
able  substitutes  for  graphical  information.  Consid¬ 
eration  has  also  be  given  to  users  with  slow  con¬ 
nections  to  the  internet. 

Using  the  Web  inside  ICASE 

The  Web  at  ICASE  is  used  to  distribute  informa¬ 
tion  to  consultants,  visitors,  and  staff  members,  as 
well  as  the  world  at  large.  Consultants  and  summer 
visitors  use  the  web  to  obtain  up-to-date  informa¬ 
tion  about  ICASE  when  they  arrive.  A  description 
of  the  local  computing  environment,  visitor  infor¬ 
mation  and  phone  numbers,  and  miscellaneous  in¬ 
formation,  such  as  paydays  and  holidays,  are  all 
available  internally  using  a  web  browser. 

The  web  provides  a  standard  access  method  that 
allows  visitors  at  ICASE  to  find  the  information 
they  need  quickly.  Visitors  come  from  a  variety  of 
computing  environments  and  use  the  information 
presented  via  the  web  to  acclimate  themselves  to 
the  computing  environment  at  ICASE. 

The  continuous  influx  of  people  at  ICASE  makes 
it  necessary  to  quickly  locate  resident  researchers. 
To  this  end,  we  are  currently  testing  an  online  form 
to  search  for  an  ICASE  staff  member  or  consultant 
and  bring  up  a  map  which  will  indicate  his/her 
location  within  the  building. 

Figure  1  shows  a  Web  page  from  the  ITRS  server 
at  ICASE  and  illustrates  the  use  of  a  search  term 
to  identify  ICASE  reports  that  contain  keywords 
of  interest  to  the  user.  The  Netscape  Web  browser 
(found  at  the  same  location  as  Mosaic)  was  used 
in  this  example.  The  report  title  is  underlined, 
denoting  that  it  is  a  Hypertext  link  to  the  report 
itself.  By  selecting  a  link,  the  user  may  download 
the  full  text  of  the  report  for  further  research. 

Usage  Statistics 

Statistics  for  the  period  May  19  through  June  20, 
1995  indicate  that  a  daily  average  of  1474  files  con¬ 
taining  a  total  of  94Mbytes  were  transferred  to  Web 
users.  This  includes  both  internal  and  external  use, 
but  it  indicates  that  the  ICASE  home  page  is  being 
heavily  used  for  its  intended  purpose. 


Fl«  Edit  Vtew  Go  Bookmorlcs  OpUans  Dtractory 


•  <?»  i: 

:  kotn 

i;  iSi  i 

Location: 

WAIS  Search  of  "compressible"  in  icase  index 


TKi*  if  k  cktrchtblemdex  Enw  keyword* 


Index  ic4<e_u>4ex  tontkun*  the  loUowmj  14  rem*  relevent  to  'eoraprM*ihle  The  fc*t  figure  tor  each  eroy 
i*  h*  relative  aoore,  the  jecond  die  numhei  of  Imei  uithe  Ken 


»  1X0  40  fc-\-ALuAfiOK  Uf-  TUfftiUlfcNCE 

HOMOGFTNEOUa  SHVlAR  Fl.Oy 

.  6C0  3g  Y;NG  IHfi  lifiASS  F;.UX  K^Crr..Y  i.N  Ttili  ANiV 

I  600  27  A.STAflLE.PEmrt^NffilliOD  WR.TTIS  COWPT.ESSDLu 
EO UATIC-NS  I-OPa j  BO-J! jPARY  CC-N&ITiaN 

I  600  51  A  PSSUDO-SK>Uh;DCOfBTIlW  i^.-nOhaHIPFORTrlE  DFlATATIONAL 

I  ao  33  A ^Fj^;ay.?-:rAncN yoR tkv: TiJKfii.ii.vijr >1  ffx coktrjruitqvjO 
HfiVN0Lns-s"iKS35  won  va'o-  eq  ii/vTinN  ci.onii 

•  ao  35  iriE 
iNCOfJITRilSSiBLE  ELOV^ 

.  '^0  31  A  CiiP.TE£lA>-l  ow'd  ATPROaCH  ViTni 
CGKtPRESiSTBLE  FLOVS 

•  400  26  A  CC;i.KAH  VAhW.  CTfittVEHEV  Mol.ll  l/A  lV.Oiv 

K>R  .VLOWg 

I  200 33 cANGNicAi -VAW ABi Po. Hiy .81  w ?f?£:rHonenR.?r?:M7y.:5:!Av  eiJlfi?. 

•  200  60  TT4E  JIAmiZniC  Em:CT  OF  COItOZRESSJEJUmMTmyt.^ 

.  100  32  EFEECT  GFrPiaSl5-L  ^SAaENTOK  il  gAAfilLir/.LM 

HIGH  -3PE5D  BC  ITJUAyt'  lA’t'EFS 

•  100  33  A  vH  r«iKrri.'K.'":i:  gb:  n  scn.vFJti; 

I  loa 37 A :-rfLT>YOPTHTc voFr!T■KK■Qwo^'£y 76/40.. dry rx/i.T;j.c -r;£f.TAyiig; 

a  100  46  PanSSllRE  IlPDATlI'K:  h>fFTHOr«?OP  TrIE  FTnvADY-STATE  PiAFDEQilATlCrK 


Figure  1:  Example  of  a  report  search  on  the  ICASE 
WEB. 


Conclusions  and  Future  Work 

The  World  Wide  Web  has  been  expanding  geomet¬ 
rically  since  ICASE  began  using  it  in  1993.  Na¬ 
tional  Public  Radio  reports  that  “about  1200  Web 
Servers  are  being  added  weekly.”  The  Web  has  be¬ 
come  an  important  method  of  sharing  information 
to  a  large  and  growing  audience  of  researchers. 

The  explosion  in  the  number  of  computers  con¬ 
nected  to  the  Web  has  created  a  large  market  for 
Web  servers  and  browsers  with  far  more  capability 
than  originally  anticipated.  An  important  exam¬ 
ple  is  the  Hot  Java  browser.  Hot  Java  adds  greater 
multi  media  capabilities.  In  order  not  to  be  “lost 
in  the  crowd,”  we  see  the  need  for  the  ICASE  home 
page  to  be  refreshed  in  terms  of  both  look  and  con¬ 
tent  on  an  annual  basis. 

Researchers  at  ICASE  have  also  been  investigat¬ 
ing  the  possibility  of  making  multi  media  ICASE 
reports  available.  This  will  he  practical  only 
through  the  use  of  the  Web.  However,  we  will 
need  to  use  the  best  available  technology  to  pro¬ 
vide  highly  compressed  image,  animation,  and  au¬ 
dio  files  to  persons  requesting  them.  A  minute  or 
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two  of  full  motion  video  with  synchronized  sound 
can  easily  exceed  lOMbytes.  I  I  ano/or 
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Infinite-Dimensional  Optimization 
Problems 

Eyal  Ariaii^ 

Introduction 

In  an  infinite-dhnensional  optimization  problem, 
the  design  variable  is  an  element  in  a  function 
space,  as  in  the  case  of  aerodynamics  shape  op¬ 
timization  problems.  In  practice,  variables  are  dis¬ 
cretized  and  the  resulting  maximum-dimensional 
problem  is  solved  niimerically.  The  size  of  the  prob¬ 
lem,  including  the  number  of  design  variables,  de¬ 
pends  on  the  level  of  refinement  of  the  discretiza¬ 
tion.  In  the  following,  we  discuss  the  analysis  and 
numerical  approach  for  solving  such  problems. 

Discretization  of  the  Design  Space 

When  solving  infinite-dimensional  optimization 
problems  numerically,  it  is  customary  to  restrict 
the  design  space  to  a  small  and  smooth  finite¬ 
dimensional  subspace,  independent  of  the  grid  re¬ 
finement.  The  major  benefit  of  this  approach  is 
in  a  substantial  reduction  in  the  size  of  the  Hes¬ 
sian,  which  greatly  reduces  the  computational  cost. 
In  some  cases  it  is  possible  to  eliminate  many  of 
the  constraints  by  a  proper  choice  of  the  design 
space  basis  functions.  On  the  other  hand,  restrict¬ 
ing  the  design  space  to  a  small  and  smooth  finite¬ 
dimensional  subspace  might  become  a  disadvan¬ 
tage  during  the  preliminary  design  phase  since  by 
doing  so,  one  determines  a  priori  the  main  features 
of  the  solution. 

The  difficulty  in  a  direct  discretization  of  the 
infinite-dimensional  design  space  is  that  the  result¬ 
ing  discrete  problem  contains  approximately  M  — 
N~~^  design  variables,  where  N  is  the  number  of 
grid  points  in  the  computational  domain,  and  d  is 
the  space  dimension.  Typically  the  Hessian,  which 
is  a  rank  M  matrix,  is  dense  and  ill-conditioned. 
This  typically  causes  slow  convergence  when  using 
otherwise  efficient  low-rank  Quasi-Newton  meth¬ 
ods  such  as  the  BFGS  (an  exception  is  the  case 
where  the  Hessian’s  spectrum  approaches  zero  con¬ 
tinuously,  i.e.,  is  rank-deficient  as  will  be  discussed 
later  in  this  article).  Our  purpose  is  to  develop  ef¬ 
ficient  numerical  methods  to  solve  such  problems 

*Eyal  Arian  received  his  Ph.D.  in  Applied  Mathematics 
from  the  Weizmann  Institute  of  Science  in  Israel  and  is  cur¬ 
rently  a  Staff  Scientist  at  ICASE.  He  joined  ICASE  in  1994. 


which  can  be  applied  without  necessarily  restrict¬ 
ing  the  design  variables  to  a  small  and  smooth  sub¬ 
space. 

Local  Fourier  Analysis  of  the  Hessian 

A  local  analysis  is  performed  to  identify  the  differ¬ 
ent  categories  of  infinite-dimensional  optimization 
problems  and  to  serve  as  a  tool  in  the  development 
of  efficient  numerical  schemes.  The  problem  can  be 
reduced  to  one  on  the  boundary  in  Fourier  space. 
High-frequency  components  of  the  errors  in  all  the 
variables  are  considered  in  the  analysis,  which  is 
a  good  approximation  apart  for  a  few  frequencies 
in  the  lower  part  of  the  spectrum.  Therefore,  the 
analysis  becomes  more  accurate  as  the  mesh  size 
decreases  since  the  upper  bound  of  the  spectrum 
on  a  given  grid,  with  mesh  size  h,  is  given  by  ^ 
(for  a  uniform  mesh  in  two  dimensions). 

The  main  goal  of  this  analysis  is  to  approximate 
the  symbol  of  the  Hessian  operator,  H  ^  in  the  high 
frequencies  near  its  minimum.  Intuitively  speak¬ 
ing,  the  different  directions  in  the  design  space 
are  spanned  by  the  basis  ,  [x  is  the  coordi¬ 
nate  along  the  boundary  of  the  designed  surface), 
and  the  eigenvalues  of  the  Hessian  are  given  by 
the  symbol  H{lo):  The  condi¬ 
tion  number  can  be  approximated  by  where 

7  is  the  highest  power  of  u  in  Qualitatively 

we  can  differentiate  between  three  classes  of  sym¬ 
bols:  monotonically  increasing  (7  >  0),  constant 
(7  =  0),  and  monotonically  decreasing  (7  <  0) 
functions  of  the  frequency  oj.  The  constant  symbol 
problem  is  well-conditioned  while  the  other  cases 
are  ill-conditioned  and  are  either  low-frequencies 
dominant  (7  <  0)  or  high-frequencies  dominant 

(7  >  0). 

Let  us  elaborate  on  the  different  cases: 

#  In  the  case  7  >  0  different  “directions”, 
have  curvatures  which  are  proportional  to 

Therefore  as  u)  varies  from  0(1)  to  the  curvature 

is  changing  dramatically,  i.e.,  from  0(1)  to  . 
Taking  the  small  disturbance  shape  optimization 
problem  as  an  example,  using  the  Euler  equations 
and  minimizing  the  pressure  distribution  on  the 
surface  to  meet  a  desired  target  distribution,  the 
leading  term  in  the  symbol  H{u))  is  proportional 
to  |cjp.  Therefore,  this  problem  is  expected  to 
be  highly  ill-conditioned  and  high-frequency  domi¬ 
nant. 

•  The  case  of  a  non-monotonic  symbol,  7  =  0, 
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is  relatively  ^^easy”  in  the  sense  that  conventional 
gradient-based  algorithms  will  have  very  good  con¬ 
vergence  properties  (the  Hessian  has  a  small  con¬ 
dition  number). 

•  The  case  of  monotonically  decreasing  symbols, 
7  <  0,  should  be  treated  with  care.  In  this  case, 
high-frequency  errors  can  enter  the  design  variables 
resulting  in  a  minor  change  in  the  cost  functional 
(or  the  gradients).  This  is  precisely  the  type  of 
problem  where  the  restriction  of  the  design  space 
to  a  small,  smooth  subspace  is  recommended.  After 
eliminating  this  difficulty,  Quasi-Newton  methods 
should  be  efficient  in  solving  this  class  of  problems, 
since  the  low-frequency  components  only  span  a 
small-dimensional  design  space. 

Note  that  the  above  analysis  should  be  per¬ 
formed  in  the  discrete  space  so  the  effects  of  a  spe¬ 
cific  numerical  scheme  on  properties  of  the  Hessian 
are  considered.  Also,  in  nonlinear  problems,  the 
non-constant  coefficients  should  be  frozen  result¬ 
ing  in  a  local  approximation,  ^(cj,xo),  where  xq 
represents  a  point  on  the  boundary. 

Efficient  Solution  using  Multigrid  Methods 

Multigrid  methods  are  particularly  attractive  when 
solving  an  infinite-dimensional  optimization  prob¬ 
lem  since  the  design  variables  are  updated  on  all 
levels  as  part  of  the  multigrid  solution  of  the  anal¬ 
ysis  and  adjoint  problems  (the  “one-shot”  method). 
The  main  difficulty  is  in  providing  a  minimiza¬ 
tion  algorithm  which  smoothes  the  design  vari¬ 
ables.  If  the  Hessian  is  ill-conditioned  and  high- 
frequency  dominant  (7  >  0),  we  expect  that  the 
gradient  descent  iteration  should  smooth  the  er¬ 
rors  in  the  design  variables.  In  elliptic  problems, 
a  high-frequency  change  of  a  quantity  which  is  de¬ 
fined  on  the  boundary  will  result  in  a  local  effect 
close  to  the  boundary;  thus  the  solution  should  be 
updated  only  in  a  small  neighborhood  close  to  the 
boundary,  resulting  in  a  low  cost  functional  evalu¬ 
ation  on  fine  levels.  Therefore,  the  solution  of  the 
optimization  problem  can  be  achieved  at  a  cost  of 
only  a  few  solutions  of  the  analysis  problems. 

The  above  ideas  were  tested  numerically  to  solve 
various  optimization  problems  including  an  opti¬ 
mal  shape  design  problem  in  two  dimensions,  us¬ 
ing  a  Poisson  equation,  where  the  design  variables 
are  the  grid  points  along  the  boundary.  The  re¬ 
sults  show  a  fast  convergence  rate  independent  of 
the  mesh  size  (See  Figure  1). 


Single  Grid  Preconditioning  Methods 

Sometimes  it  is  desirable  to  have  a  method  which 
can  be  applied  when  the  analysis  code  is  using  one 
level  of  discretization.  In  that  case,  local  Fourier 
analysis  can  be  used  to  compute  a  preconditioned 
direction,  WuF^  which  is  a  good  approximation 
for  the  Newton  direction  in  the  high-frequencies. 
The  preconditioner,  V,  is  computed  such  that  its 
symbol  approximates  the  inverse  of  the  Hessian  in 
the  high-frequency  regime:  V{uj)  ^  1). 

In  practice,  the  preconditioned  direction  should  be 
projected  onto  the  finite-dimensional  design  space. 

The  proposed  algorithm  uses  two  directions  in 
each  minimization  step  as  is  done  in  the  trust  re¬ 
gion  double  dogleg  algorithm.  The  two  directions 
are  given  by  the  preconditioned  direction,  WuF\ 
and  a  Quasi-Newton  direction  which  accounts  for 
the  lowest  frequencies.  This  method  is  currently 
under  research. 
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Figure  1.  Convergence  Rates  for  an  Optimal- 
Shape  Design  Problem  Using  Multigrid  Methods. 
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A  Model  for  Velocity  Gradients 
in  Homogeneous  Turbulence 

Sharath  S.  Girimajt* 

Charles  G.  Speziale^ 


the  material  presented  here,  please  see  Girimaji 
and  Speziale  [3].  The  modified  restricted  Euler 
equation  for  the  total  velocity  gradients  is  given 
by: 


Introduction  and  Modeling 

The  phenomenon  of  turbulence,  it  is  generally 
agreed,  is  completely  described  by  the  Navier- 
Stokes  equations.  However,  the  Navier-Stokes 
equations  have  no  known  general  analytical  solu¬ 
tion,  and,  even  numerical  solutions  are  possible 
for  only  the  simplest  of  flow  geometries  at  low 
Reynolds  numbers.  The  objective  of  turbulence 
modeling  is  to  devise  an  alternative  set  of  tractable 
equations  that  mimic  certain  aspects  of  the  turbu¬ 
lence  phenomenon.  Most  models  are  phenomeno¬ 
logical  in  nature  and  attempt  to  make  quantitative 
predictions  of  turbulent  flow  statistics.  There  is 
another  genre  of  models  whose  objective  is  to  elu¬ 
cidate  some  physics  of  the  turbulence  phenomena 
without  necessarily  being  quantitatively  accurate. 
One  such  model  is  the  restricted  Euler  equation 
(Vieillefosse  [1],  Cantwell  [2])  given  by 

^  +  AikAkj  -  lA,nnAn,nSiJ  =  0.  (1) 

at  3 

where  Aij  is  the  velocity  gradient  tensor  and  Sij 
is  the  Kronecker  delta  function.  This  model  equa¬ 
tion  reproduces  some  important  turbulent  small- 
scale  features  frequently  observed  in  direct  numer¬ 
ical  simulations  (DNS)  of  isotropic  and  homoge¬ 
neous  shear  flows.  However,  when  the  restricted 
Euler  equation  is  compared  against  DNS  databases 
for  more  complex  turbulent  flows,  the  agreement 
is  not  as  good,  especially  in  regions  of  moderate 
to  low  dissipation.  This  shortcoming  may  be  due 
to  the  fact  that  it  violates  the  balance  of  mean 
momentum  for  virtually  all  homogeneous  turbulent 
flows  with  constant  mean  velocity  gradients  except 
isotropic  and  homogeneously  sheared  turbulence. 

Recently,  we  have  developed  a  modified  version 
of  the  restricted  Euler  equation  which  satisfies  the 
mean  momentum  equation  for  general  turbulent 
flows  and,  hence,  is  better  capable  of  simulating 
more  complex  test  cases.  For  a  detailed  report  of 
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dA 


dt 


T  AikAf;j 


r  AjYifi  Ajiiji  Sij 


(2) 


dA-  —  —  1—  — 

-p  AikA^j  —  -ijAmnAnm^ij- 
dt  o 


A  comparison  of  equations  (1)  and  (2)  reveals  the 
difference  between  the  original  and  the  modified  re¬ 
stricted  Euler  equation.  In  the  modified  restricted 
Euler  equation  model,  the  mean  velocity  gradient 
terms  on  the  right-hand-side  of  equation  (2)  sur¬ 
vive  since  the  anisotropic  mean  pressure  Hessian  is 
not  neglected  —  a  feature  that  maintains  consis¬ 
tency  with  the  balance  of  mean  momentum.  The 
original  restricted  Euler  equation  is  recovered  for 
isotropic  turbulence  and  homogeneous  shear  flow. 

Results  and  Discussion 

We  solve  the  modified  restricted  Euler  equation 
numerically  for  isotropic  and  anisotropic  homoge¬ 
neous  turbulent  flows.  An  ensemble  of  4096  fluid 
particles  is  considered  and  the  velocity  gradient 
for  each  particle  is  obtained  numerically  by  inte¬ 
grating  the  modified  restricted  Euler  equation  us¬ 
ing  a  predictor-corrector  Runge-Kutta  scheme.  We 
present  the  behavior  of  the  the  anisotropic  struc¬ 
ture  tensor  function  defined  as 

_  1 

Fij  —  hi^ibfi^j  -^Sij ,  (3) 

In  the  above  equation  bij  =  Aij J \/ AmnAmn-  In 
isotropic  turbulence,  the  components  of  these  ten¬ 
sors  are  all  zero.  The  magnitudes  and  signs  of 
the  non-zero  components  should  therefore  convey  a 
good  description  of  the  anisotropic  structure  of  the 
fluctuating  velocity  gradient  tensor.  We  also  exam¬ 
ine  some  intrinsic  properties  of  the  velocity  gradi¬ 
ents.  The  results  presented  here  are  for  the  case  of 
a  rapid  solid  body  rotation:  Aij  =  u)Sii8j2—a)Si2Sj\. 

The  time  evolution  of  the  diagonal  components 
of  the  anisotropic  structure  tensor  function  Fij  is 
given  in  Figure  1.  Rapid  solid  body  rotation  of 
an  initially  isotropic  velocity  field  does  not  appear 
to  affect  the  isotropy  of  the  structure  tensor  func¬ 
tion  immediately.  AU  three  components  stay  close 
to  zero  initially  for  a  short  duration  of  time.  This 
behavior  is  consistent  with  the  observations  from 
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Figure  1:  Time  evolution  of  Fij  (new  model). 

the  DNS  of  an  initially  isotropic  velocity  field  sub¬ 
jected  to  a  rapid  solid  body  rotation.  The  DNS 
results  indicate  that  rapid  rotation  produces  iner¬ 
tial  waves  that  scramble  the  turbulence  and  sur¬ 
vive  for  several  eddy  turnover  times  leaving  the  ve¬ 
locity  field  with  approximately  isotropic  one-point 
statistics.  Once  these  inertial  waves  are  damped, 
anisotropy  due  to  mean  rotation  can  set  in.  At 
later  times,  the  axial  component,  F33,  starts  to  de¬ 
crease  (in  comparison  with  Fu  and  F22)  gradually 
and  monotonically,  a  feature  that  is  consistent  with 
the  Taylor- Proudman  theorem  which  calls  for  the 
axial  velocity  gradient  to  vanish  at  infinite  rota¬ 
tion  rate.  The  original  restricted  Euler  equation, 
on  the  other  hand,  predicts  that  the  axial  compo¬ 
nent  of  the  structure  function  tensor,  F33,  increases 
with  time,  completely  inconsistent  with  the  Taylor- 
Proudman  theorem  (Figure  2). 


Figure  2:  Time  evolution  of  Fij  (old  model). 

Intrinsic  geometry  of  bij. 

Examination  of  DNS  databases  reveal  that  the  in¬ 
trinsic  geometry  of  the  fluctuating  velocity  gradi¬ 
ents  in  anisotropic  homogeneous  turbulence  is  sim¬ 


ilar  to  that  in  isotropic  turbulence.  The  eigenvec¬ 
tors  of  the  strain-rate  tensor  are  ai,  a2  and  as, 
where  the  eigenvalues  are  arranged  in  descending 
order.  Due  to  incompressibiUty,  which  renders  the 
divergence  of  the  strain  rate  rate  zero,  the  largest 
eigenvalue  is  constrained  to  be  positive  and  the 
smallest  one  negative.  The  modified  restricted  Eu¬ 
ler  equation,  consistent  with  the  original  model  and 
DNS  findings,  predicts  that  the  vorticity  vector 
is  best  oriented  with  the  intermediate  eigenvector 
(a2).  The  evolution  of  the  averaged  eigenvalues 
is  given  in  Figure  3.  The  averaged  intermediate 
eigenvalue  is  indeed  positive,  again  consistent  with 
DNS  findings.  Furthermore,  these  intrinsic  prop¬ 
erties  are  also  preserved  for  both  the  homogeneous 
shear  flow  and  plane  strain  flow  cases  which  are  not 
shown. 

These  results  lead  us  to  believe  that  the  new 
modified  restricted  Euler  equation  can  serve  as  a 
valuable  dynamical  model  for  the  analysis  of  the 
geometry  of  the  velocity  gradient  tensor  in  general 
homogeneous  turbulent  flows. 


Figure  3:  Time  evolution  of  averaged  eigenvalues. 
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Numerical  Simulations  of 
Shock- Vortex  Interactions 

Gordon  Erlebacher* 

Thomas  L.  Jackson^ 

In  the  late  1950’s,  Ribner  began  his  study  of 
shock/turbulence  interactions.  To  this  end,  he  con¬ 
sidered  a  turbulent  flow  as  a  collection  of  linear 
waves  which  could  be  superimposed  to  form  a  re¬ 
sulting  turbulent  field.  By  studying  the  properties 
of  a  single  wave  interacting  with  a  shock,  he  was 
able  to  deduce  the  effects  of  a  shock  on  an  upstream 
turbulent  flow.  Since  then,  there  have  been  sev¬ 
eral  studies  seeking  to  extend  Ribner’s  work.  Thus, 
upstream  acoustic  and  entropy  waves  were  consid¬ 
ered,  and  experimental  work  was  initiated  both  to 
corroborate  the  theory,  and  extend  it.  Ten  years 
ago,  the  first  simulations  were  performed  by  Hus- 
saini  and  his  collaborators  [1,2].  First  with  finite- 
difference  methods,  then  spectral  methods,  and 
more  recently  back  to  high-order  finite-difference 
methods,  their  studies  mostly  addressed  the  pa¬ 
rameter  range  under  which  linear  theory  was  ex¬ 
pected  to  hold.  Several  nonlinear  simulations  of 
upstream  entropy  spots,  vortex  pairs  and  vortex 
streaks  interacting  with  a  shock  were  performed  as 
well.  However,  downstream  boundary-conditions 
were  found  to  affect  the  results  and  thus,  long¬ 
time  integrations  could  not  be  performed.  Further¬ 
more,  the  quantitative  study  of  nonlinear  effects 
were  never  performed. 

A  related  problem,  which  can  be  viewed  as  an 
important  element  of  the  more  complex  process  of 
shock-turbulence  interaction,  is  the  interaction  of  a 
single  vortex  and  a  shock  wave.  This  process  is  also 
of  interest  in  aeroacoustics  owing  to  the  production 
of  noise  generated  by  shocks  present  in  high-speed 
aircraft.  Thus,  in  the  second  part  of  this  work,  the 
interaction  of  a  single  vortex  and  a  shock  wave  is 
also  investigated  and  compared  to  linear  results. 
Results  are  presented  for  various  vortex  strengths 
and  at  various  Mach  numbers. 

Although  at  the  time,  the  primary  motivation 
behind  these  studies  was  to  understand  the  various 
types  of  amplification  mechanisms  of  turbulence 
passing  through  a  shock,  recent  work  has  focused 
on  the  acoustic  properties  of  the  sound  generated 
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when  a  vortex  crosses  a  shock.  Indeed,  a  single 
vortex  produces  an  acoustic  wave  which  propagates 
away  from  the  vortex  core  at  the  speed  of  sound, 
and  is  responsible  for  the  far-field  noise  measured 
by  experimentalists. 

We  consider  an  ideal  gas  governed  by  the  two- 
dimensional  Euler  equations  in  nonconservation 
form,  supplemented  with  the  equation  of  state  for 
an  ideal  gas.  The  Euler  equations  are  solved  us¬ 
ing  a  shock-fitted  algorithm.  To  this  end,  two 
separate  computational  domains  are  required:  one 
for  the  upstream  region,  and  one  for  the  down¬ 
stream.  Rather  than  construct  a  computational 
domain  which  abuts  the  downstream  domain,  we 
utilize  the  fact  that  the  upstream  (supersonic  rela¬ 
tive  to  the  shock)  flow  is  not  affected  by  the  shock 
motion.  Therefore,  all  that  is  required  to  prop¬ 
erly  determine  the  downstream  solution  (including 
the  shock  position)  is  the  upstream  solution  at  the 
shock.  In  this  way,  the  upstream  solution  can  be 
calculated  independently  of  the  downstream  solu¬ 
tion.  An  evolution  equation  for  the  shock  is  de¬ 
rived  by  applying  the  Rankine-Hugoniot  conditions 
to  the  downstream  characteristic  pointing  into  the 
shock,  at  the  shock  itself.  Thus,  both  shock  posi¬ 
tion  and  shape  evolve  in  time.  All  spatial  deriva¬ 
tives  are  computed  using  a  compact  order  dis¬ 
cretization  scheme,  with  special  treatment  at  the 
boundaries  to  maintain  6^'’-  order  accuracy.  Time 
advancement  is  based  on  a  low-storage,  5  stage, 
4*'*  order  accurate  Runga-Kutta  scheme.  Instead 
of  implementing  a  characteristic  boundary  condi¬ 
tion  downstream,  we  choose  instead  to  construct  a 
buffer  domain.  In  this  buffer,  the  Euler  equations 
are  modified  to  add  a  convective  term  to  the  time 
derivative  whose  role  is  to  curve  the  characteristics 
in  such  a  way  that  the  characteristic  velocities  are 
all  supersonic  at  the  exit  plane.  Complete  details 
are  given  in  [2]. 

Although  many  studies  have  been  made  of 
shock/vortex  interactions,  there  have  been  few 
quantitative  studies  of  the  nonlinear  processes.  It 
is  known  that  for  weak  shocks,  as  the  strength  of 
the  vortex  increases,  the  shock  develops  a  triple 
point,  which  has  also  been  observed  experimen¬ 
tally.  However,  there  are  instances  where  a  triple 
point  does  not  form  and  nonlinear  processes  are  ac¬ 
tive  nonetheless.  To  address  this  issue,  we  place  a 


8 


vortex  with  azimuthal  velocity  profile 


_ - _  ^  ~0.5(r/a)^ 

27ra^e-o^^ 


in  the  upstream  domain.  The  pressure,  density 
and  temperature  satisfy  the  isentropic  relation  p  = 

7-1 

p^‘  =  T  ^  .  The  analytical  form  of  this  vortex  is 
an  exact  solution  to  the  nonlinear  Euler  equations. 
We  choose  a  core  radius  a  =  0.3,  and  study  the 
effect  of  circulation  F  and  shock  Mach  number  M 
on  the  flow  downstream  of  the  shock.  The  radial 
exponential  decay  of  all  variables  insures  that  spu¬ 
rious  reflections  will  not  occur  from  the  freestream 
boundaries. 

To  identify  nonlinear  effects,  we  computed  the 
RMS  kinetic  energy,  pressure  energy,  vorticity  and 
dilatation  in  the  downstream  domain,  along  with 
the  dominant  terms  which  contribute  to  their  rate 
of  change  with  time.  All  variables  were  normal¬ 
ized  with  r  to  insure  that  nonlinear  effects  mani¬ 
fest  themselves  as  a  departure  from  the  time  his¬ 
tories  at  infinitesimal  F.  We  found  that  the  dilata- 
tional  energy  quickly  rises  from  zero  to  reach  some 
maximum  value  shortly  after  the  vortex  crosses  the 
shock.  For  moderate  F,  this  rise  is  followed  first 
by  a  short  period  of  fast  decay,  which  corresponds 
to  the  formation  of  the  acoustic  pulse  and  its  ini¬ 
tial  separation  from  the  vortex  core.  The  ensu¬ 
ing  slower  decay  is  probably  due  to  energy  loss 
as  the  pulse  propagates  outward  from  the  vortex 
center.  As  F  Increases,  the  pressure  and  dilata¬ 
tion  “saturation”  levels  scale  approximately  as  F^. 
The  scaling  is  not  exact  because  a  contribution  of 
pressure  is  present  due  to  the  vortical  field.  It 
is  clear  that  the  F^  scaling  becomes  more  appro¬ 
priate  as  the  vorticity  becomes  stronger,  and  the 
downstream  contribution  of  upstream  pressure  be¬ 
gins  to  dominate  that  of  the  pressure  due  to  the 
upstream  vortical  modes.  Note  that  a  theory  ca¬ 
pable  of  making  this  prediction  must  be  at  least 
taken  one  step  beyond  linear.  We  have  performed 
a  high  Mach  number  asymptotic  expansion  of  the 
linearized  Rankine-Hugoniot  conditions,  followed 
by  a  second  expansion  in  the  vortex  circulation. 
The  analysis  reveal  that  nonlinear  effects  domi¬ 
nate  the  flow  once  FM  >  constant.  Moreover,  the 
quadratic  dependence  of  the  downstream  perturba¬ 
tion  pressure  on  the  shock  Mach  number  is  not  the 
result  of  a  finite-pressure  inside  the  upstream  vor¬ 
tex,  but  rather  is  due  to  the  nonzero  density  inside 


the  vortex,  due  to  imposed  isentropic  conditions. 
Finally,  as  the  circulation  is  increased  to  F  =  1,  we 
note  that  the  dilatational  energy  begins  to  increase 
in  time.  We  hypothesize  that  this  is  due  to  the 
steepening  of  the  acoustic  pulse  as  it  radiates  away 
from  the  vortex  core.  This  is  further  evidenced 
by  the  two-dimensional  plots  shown  in  Figure  1. 
Another  effect  of  nonlinearity  is  the  destruction  of 
the  quadrupole  structure  of  the  acoustic  pulse.  As 
seen  from  the  low  F  contours,  the  pressure  field 
(and  the  dilatation)  have  a  quadrupole  structure. 
The  two  negative  lobes  have  merge  as  the  circula¬ 
tion  increases,  and  the  symmetry  with  respect  to 
the  y  =  0  axis  has  disappeared.  More  figures  and 
details  are  available  in  [2]. 


X 


X 


Figure  1:  Contour  lines  of  pressure.  Mach  8  shock, 
F  =  0.03  (top)  and  F  =  1.0  (bottom).  Negative 
values  are  dotted,  positive  values  are  sobd. 
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June  1. 

McHugh,  P.  R.,  Idaho  National  Engineering  Laboratory,  “Inves¬ 
tigation  of  Newton-Krylov  Algorithms  for  Solving  hicom- 
pressible  and  Low  Mach  Number  Compressible  Flow  Using 
Finite  Volume  Discretization,”  June  1. 

Barber,  Thomas,  United  Technologies  Research  Center,  “Turbu¬ 
lence  Modeling  for  Aerospace  Commercial  Application,” 
June  2. 

Clriueh,  Tzi-cker,  State  University  of  New  York  at  Stony  Brook, 
“Data  Compression:  A  Computer  Systems  Perspective,” 
June  6. 

Hanebutte,  Ulf,  hitelligent  Transportation  Systems,  Argomie 
National  Laboratory,  “Parallel  Traffic  Simulations  —  An 
Overview,”  Jmie  7. 


Yoke,  Peter  R.,  University  of  Surrey,  England,  “Large-Eddy 
Simulation  of  Boundary  Layer  Bypass  Transition  Mecha¬ 
nisms,”  June  8. 

Crimiiiale,  William,  University  of  Washington,  “Disturbance 
Transient  Dynamics  in  Shear  Flows,”  June  15. 

Yeung,  P.K.,  Georgia  Institute  of  Technology,  “Parallel  Compu¬ 
tations  of  Turbulent  Mixing  and  Dispersion,”  June  19. 

Reuther,  James,  RIACS,  “Practical  Aspects  of  Aerody¬ 
namic  Shape  Optimization  Using  Adjoint  Based  Sensitivity 
Derivatives,”  June  20. 

Nordstrom,  Jan,  The  Aeronautical  Research  histituteof  Sweden, 
“The  Use  of  Characteristic  Boundary  Conditions  for  the 
Navier-Stokes  Equations,”  June  26. 

Rao,  Jagaimatha,  University  of  Houston,  “New  Models  for  Opti¬ 
mal  Truss  Topology  in  Limit  Design  Based  on  Unified  Elas¬ 
tic/Plastic  Analysis,”  June  28. 

McMichael,  James,  Air  Force  Office  of  Scientific  Research, 
“AFOSR  Research  in  Turbulence  and  hiternal  Flows,” 
June  29. 

ICASE  REPORTS 

April  1,  1995  -  June  30,  1995 

Dando,  Andrew,  and  Philip  Hall:  Wavenumber  selection  for 
small-wavelength  Goriler  vortices  in  curved  channel  flows. 
ICASE  Report  No.  95-21,  April  6,  1995,  28  pages.  Submitted 
to  Studies  in  Applied  Mathematics. 


Ristorcelli,  J.  R.:  A  pseudo-sound  constitutive  relationship  for 
the  dilatational  covariances  in  compressible  turbulence:  An  ana¬ 
lytical  theory.  ICASE  Report  No.  95-22,  April  4,  1995,  50  pages. 
Submitted  to  the  Journal  of  Fluid  Mechamcs. 


Shortis,  Trudi  A.,  and  Phihp  HaU:  On  the  nonlinear  stability  of 
the  unsteady,  viscous  flow  of  an  incompressible  fluid  in  a  curved 
pipe.  ICASE  Report  No.  95-23,  April  6,  1995,  23  pages.  Submit¬ 
ted  to  the  Journal  of  Fluid  Mechanics. 


Sun,  Xian-He,  and  Jianping  Zhu:  Performance  prediction:  A  case 
study  using  a  multi-ring  KSR-1  machine.  ICASE  Report  No.  95- 
24,  April  6,  1995,  21  pages.  To  appear  in  the  Proceedings  of  the 
28th  Hawaii  International  Conference  on  System  Sciences,  Jan. 
’95. 


Haras,  Zigo,  and  Shlomo  Ta’asan:  The  large  discretization  step 
method  for  time- dependent  partial  differential  equations.  ICASE 
Report  No.  95-25,  April  7,  1995,  53  pages.  To  be  submitted  to 
the  Journal  of  Computational  Physics. 


Mavriphs,  D.  J.:  Unstructured  mesh  generation  and  adaptivity. 
ICASE  Report  No.  95-26,  April  13,  1995,  47  pages.  Lecture  notes 
for  26th  CFD  Lecture  Series  of  von  Karman  Institute;  AGARD 
Pubhcation. 
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Mavriplis,  D.  J.:  Multigrid  techniques  for  unstructured  meshes. 
ICASE  Report  No.  95-27,  April  13,  1995,  61  pages.  Lecture  notes 
for  26th  CFD  Lecture  Series  of  von  Karman  Institute;  AGARD 
Publication. 


Venkatakrisiman,  V,:  hnplicit  schemes  and  parallel  computing 
in  unstructured  grid  CFD.  ICASE  Report  No.  95-28,  April  13, 
1995,  62  pages.  Lecture  notes  for  26th  CFD  Lecture  Series  of  von 
Karman  Institute. 


Hall,  Philip,  and  Deinetrios  T.  Papageorgiou:  On  the  modula- 
tional  instability  of  large  amplitude  waves  in  supersonic  bound¬ 
ary  layers.  ICASE  Report  No.  95-29,  April  13,  1995,  44  pages. 
Submitted  to  SIAM. 


Mavriplis,  D,  J.,  and  V.  Venkatakrisiman:  A  3d  agglomeration 
■multigrid  solver  for  the  reynolds-averaged  Navier-Siokes  equa¬ 
tions  on  unstructured  meshes.  ICASE  Report  No.  95-30,  .4pril 
13,  1995,  29  pages.  Submitted  to  the  International  Journal  for 
Numerical  Methods  in  Fluids. 


Crockett,  Thomas  W.:  Parallel  rendering.  ICASE  Report  No. 
95-31,  May  10,  1995,  52  pages.  To  appear  in  “The  Encyclopedia 
of  Computer  Science  and  Technology.” 


Ta’asan:  Pseudo-time  methods  for  constrained  optimization  prob¬ 
lems  governed  by  PDF.  ICASE  Report  No.  95-32,  May  18,  1995, 
15  pages.  To  be  submitted  to  SIAM  Journal  on  Optimization. 


Stoica,  Ion,  Hussein  Abdel- Wahab,  and  Alex  Pothen:  A  microe¬ 
conomic  scheduler  for  parallel  computers.  ICASE  Report  No. 
95-33,  June  12,  1995,  23  pages.  To  appear  in  Proceedings  of  the 
Workshop  on  Job  Scheduling  Strategies  for  Parallel  Processing. 
9th  International  Parallel  Processing  Symposium. 


Morano,  E.,  D.  J.  Mavriplis,  and  V.  Venkatakrisiman:  Coarsening 
strategies  for  unstructured  multigrid  techniques  with  application 
to  anisotropic  problems.  ICASE  Report  No,  95-34,  May  31,  1995, 
19  pages.  To  be  submitted  to  SIAM  Journal  on  Scientific  Com¬ 
puting. 


Fahringer,  Thomas,  Matthew  Haines,  and  Piyush  Mehrotra:  On 
the  utility  of  threads  for  data  parallel  programming.  ICASE  Re¬ 
port  No.  95-35,  May  10,  1995,  14  pages.  To  appear  in  Interna¬ 
tional  Conference  on  Supercomputing  ’95. 


Haines,  Matthew,  Piyush  Mehrotra,  and  David  Cronk:  Ropes: 
Support  for  collective  operations  among  distributed  threads. 
ICASE  Report  No.  95-36,  May  10,  1995,  19  pages.  Submitted 
to  Supercomputing  ’95. 


Carpenter,  Mark  H.,  and  David  Gottlieb:  Spectral  methods  on 
arbitrary  grids.  ICASE  Report  No.  95-37,  May  10,  1995,  29 
pages.  To  be  submitted  to  the  Journal  of  Computational  Physics. 

D’Ambrosio,  Domenic,  and  Roberto  Marsilio:  A  n-umerical 
method  for  solving  the  three-dimensional  parabolized  Navier- 
Siokes  equations.  ICASE  Report  No.  95-38,  May  10,  1995,  24 
pages.  To  be  submitted  to  the  Journal  of  Computational  Physics. 


Bauer,  F.,  L.  Maestrello,  and  L.  Ting:  Acoustic  field  in  unsteady 
moving  media.  ICASE  Report  No.  95-39,  May  10,  1995,  47  pages. 
Submitted  to  JASA. 


Fiterman,  A.,  E.  Turkel,  and  V.  Vatsa:  Pressure  updating  methods 
for  the  steady-state  fluid  equations.  ICASE  Report  No.  95-40, 
May  10,  1995,  15  pages.  To  appear  in  the  Proceedings  of  the 
AIAA  CFD  Conference. 


Bayliss,  A.,  L.  Maestrello,  J.  L.  McGreevy,  and  C.  C.  Fenno: 
Response  of  a  multi-panel  assembly  to  noise  from  a  jet  in  forward 
motion.  ICASE  Report  No.  95-41,  May  18,  1995,  29  pages.  To 
appear  at  First  CEAS/AIAA  Aeroacoustics  Conference. 

Banks,  H.  T.,  R.  C,  Smith,  D.  E..  Brown,  R.  J.  Silcox,  V.  L. 
Metcalf:  Experimental  confirmation  of  a  PDE-based  approach  to 
design  of  feedback  controls.  ICASE  Report  No.  95-42,  May  19, 
1995,  42  pages.  Submitted  to  SIAM  Journal  of  Control  and  Op¬ 
timization. 


Zhou,  Ye:  A  phenomenological  treatment  of  rotating  turbulence. 
ICASE  Report  No.  95-43,  May  18,  1995,  13  pages.  To  appear  in 
Physics  of  Fluids, 


Girimaji,  Sharath  S.,  and  Ye  Zhou:  Analysis  and  modeling  of 
subgrid  scalar  mixing  using  numerical  data.  ICASE  Report 
No.  95-44,  May  31,  1995,  31  pages.  To  be  submitted  to  Physics 
of  Fluids. 


Burns,  John  A.,  and  Belinda  B.  King:  Representation  of  feedback 
operators  for  hyperbolic  systems.  ICASE  Report  No.  95-45,  May 
31 ,  1995,  13  pages.  To  appear  in  the  “Conference  on  Computation 
and  Control,”  J.  Lund,  Ed.,  Birkhauser,  1995. 

Middleton,  David:  An  executable  specification  for  the  message 
processor  in  a  simple  combining  network.  ICASE  Report  No.  95- 
47,  June  12,  1995,  36  pages. 

Tidriri,  M.  D.:  Krylov  methods  for  compressible  flows.  ICASE 
Report  No.  95-48,  June  30,  1995,  50  pages.  To  be  submitted  to 
Computer  Methods  in  Applied  Mechanics  and  Engineering. 

Bokhari,  Shahid  H.:  The  LINUX  operating  system:  An  introduc¬ 
tion.  ICASE  Report  No.  95-49,  June  30,  1995,  20  pages.  To 
appear  in  IEEE  “Computers” . 


Banks,  David  C.,  and  Michael  Kelley:  Tracking  a  turbulent  spot  in 
an  immersive  environment.  ICASE  Report  No.  95-50,  Jmie  30, 
1995,  12  pages.  Submitted  to  the  1995  Symposium  on  Interactive 
3D  Computer  Grapliics. 


Banks,  David  C.,  and  Bart  A.  Singer:  A  predictor- corrector  tech¬ 
nique  for  visualizing  unsteady  flows.  ICASE  Report  No.  95-51, 
Jmie  30,  1995,  30  pages.  Submitted  to  IEEE  Transactions  on 
Visualization  and  Computer  Graphics. 

Macrossan,  M.  N.:  Some  developments  in  the  equilibrium  particle 
simulation  method  for  the  direct  simulation  of  compressible  flows. 
ICASE  Interim  Report  No.  27,  June  27,  1995,  33  pages. 
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Upcoming  Symposium 

Graphics  and  visualization  is  an  important  com¬ 
ponent  of  the  computer  science  research  program 
at  ICASE.  To  update  our  knowledge  in  this  area, 
ICASE/LaRC  will  hold  a  symposium  this  year. 

•  Symposium  on  Visualizing  Time-Varying  Data 
—  September  18-19,  1995 

Further  information  may  be  found  on  the  World 
Wide  Web  at 

http:// www.ic.ase.edu  /  workshops/ vtvd. 

Publication  Requests 

If  you  are  interested  in  receiving  specific  ICASE 
reports,  either  write  or  email  your  requests  to: 

Shannon  Keeter 

ICASE 

Mail  Stop  132C  ■ 

NASA  Langley  Research  Center 
Hampton,  VA  2^81-0.001 
pubs@icase.edu 

Please  allow  four  weeks  for  delivery. 


ICASE 

Mail  Stop  1 32C 

NASA  Langley  Research  Center 
Hampton,  VA  23665 


The  ICASE  Research  Quarterly  is  published  by 
the  Institute  for  Computer  Applications  in  Science 
and  Engineering,  The  Institute  is  an  affirmative  ac¬ 
tion/equal  opportunity  employer. 

Editor 

Gordon  Erlebacher 
Editorial  Assistants 

Emily  Todd  and  Shelly  Millen 


Further  Information 

Further  information  about  ICASE  or  the  contents 
of  this  newsletter  can  be  obtained  by  contacting: 

Emily  Todd 

ICASE 

Mail  Stop  132C 

NASA  Langley  Research  Center 
Hampton,  VA  23681-0001 
Fax:  (804)  864-6134 
info@icase.edu 
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